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Abstract 

We consider General Relativity with matter, radiation and a min¬ 
imally coupled dark energy defined by an equation of state w. Using 
dynamical system method, we find the equilibrium points of such a 
theory assuming an expanding Universe and a positive dark energy 
density. Two of these points correspond to classical radiation and 
matter dominated epochs for the Universe. For the other points, dark 
energy mimics matter, radiation or accelerates Universe expansion. We 
then look for possible sequences of epochs describing a Universe start¬ 
ing with some radiation dominated epoch(s) (mimicked or not by dark 
energy), then matter dominated epoch(s) (mimicked or not by dark en¬ 
ergy) and ending with an accelerated expansion. We find ten sequences 
able to follow this Universe history without singular behaviour of w at 
some saddle points. Most of them are new in dark energy literature. To 
get more than these ten sequences, w has to be singular at some specific 
saddle equilibrium points. This is an unusual mathematical property 
of the equation of state in dark energy literature, whose physical conse¬ 
quences tend to be discarded by observations. This thus distinguishes 
the ten above sequences from an infinity of ways to describe Universe 
expansion. 


1 Introduction 

An important problem of today cosmology is to choose between an infinite 
number of theories which ones are able to explain observations. Most of 
time, this is done by comparing a theory or a whole family of theories with 
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some observations as supernovae, CMB, weak lensing and so on. Another 
basic, and at first glance obvious question to validate a cosmological theory, 
is also to check if it is able to reproduce a Universe history starting by a 
radiation dominated epoch, followed by a matter dominated epoch and then 
an accelerated expansion epoch. For instance it has been claimed that this is 
not possible for f(R) theories in the metric formalism(T] but if they respect 
some special conditions, a contested result in [2]. In this paper, we exam¬ 
ine how General Relativity with matter, radiation and a minimally coupled 
dark energy defined by an equation of state w is able to reproduce the above 
Universe history whose ACDM description is one way among many others. 
To reach this aim, we apply dynamical system method to the equations of 
General Relativity, rewritten with a set of normalised variables. We get the 
equilibrium points of this equations system and give some general (since 
we will keep w unspecified) considerations about their stability. Two of 
these points represent usual radiation and matter dominated epochs. For 
the other points, dark energy mimics matter or radiation but for one which 
corresponds to an accelerated expansion of the Universe. Then we try to de¬ 
termine all the possible trajectories in the phase space describing sequences 
of epochs starting by some dominated radiation epochs (mimicked or not 
by dark energy), then some matter dominated epochs (mimicked or not by 
dark energy) and ending with an accelerated expansion epoch. We find ten 
of these sequences when w is not singular at some saddle points. Most of 
them are new in the dark energy literature like the ones when dark energy 
mimics radiation. To get more than these ten sequences, a necessary but not 
sufficient condition is that w be singular at some saddle points. Physically, 
since a trajectory in the phase space never pass through a saddle point but 
close to it, it means that w should have pronounced bumps or/and dips, a 
property of w that tends to be discarded by observations jT3] [T4| . 

These results may be seen as a classification of dark energy theories in 
presence of matter and radiation, based on dynamical system. There are 
already several dark energy classifications in the literature. In [3], dark en¬ 
ergy theories are classified in a plane of observables that corresponds to the 
parameterisation of a variable equation of state, w(a) = wo + w a ( 1 — a), 
where a is the scale factor of the universe. A similar study has been made 
in [5j. Another way of classifying dark energy models is using statefinder 
parameters [9]. In |3], the class of all dark energy models is considered as 
a multiverse and classified using the notion of Banach space. The classifi¬ 
cation we get in this paper is different from these ones. It is based on the 
notion of matter, radiation and dark energy equilibrium points but it does 
not distinguish between a quintessence or a ghost model for instance. 

The plan of this paper is as follows. In the second section, we look for all 
possible equilibrium points of the dynamical system defined by equations of 
General Relativity with matter, radiation and a minimally coupled dark en¬ 
ergy. We give some general considerations about their stability. In the third 
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section, we look for some possible sequences of epochs describing Universe 
first dominated by radiation, then matter and then dark energy accelerating 
expansion. We find 10 of these sequences when w is not singular at some 
saddle equilibrium points. We discuss our results in the last section. We 
explain why w should be singular at some saddle equilibrium points to re¬ 
produce an infinity of sequences and not only the 10 of the third section. 
This property of w could be in tension with the data, thus distinguishing 
the 10 sequences we found in the third section among all the possible ones. 


2 Field equations 


Field equations of General Relativity with a dark energy in a flat FLRW 
Universe write 


dH 

dt 


k 

= — {p m + Pr + P) 

k 

- [-3 Pm - 4 Pr - 3(tt> + 1 )p\ 


p m is the matter density. From conservation of the stress-energy tensor for 
p m we derive that p rn = p m o(l + z) 3 , with z the redshift and p m o the matter 
density today. p r is the radiation density. From conservation of the stress- 
energy tensor for p r we derive that p r = p r o(l + z) A with p r o the radiation 
density today, p is the dark energy density with equation of state p/p = w, 
p being the dark energy pressure. Conservation of the stress-energy tensor 
for p writes 


Then we define 


= —3 H{w + 1 )p 

(1) 

k Pm 

Vl ~ 3 IP 

(2) 

k p r 
m ~ 3 IP 

(3) 

k p 
m ~ 3 IP 

(4) 


With the constraint equation 


l = yl + y 2 + y 3 

We are going to make two assumptions: 

1. H > 0, i.e. Universe is expanding in agreement with the data. 

2. p > 0, i.e. dark energy density is positive. Negative energy density 
generally belongs to the realm of quantum cosmology that we will not 
consider in this work. 
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Since p m > 0, p r > 0 and we assume p > 0, the above equation implies that 
0 < Hi < 1. Differentiating 2/1 and y 2 and considering field equations and 
energy density conservation of p m , p r and p, we are left with the following 
dynamical system 



A prime stands for a derivative with respect to N = In a, with a the scale 
factor. Since we assume that Universe is expanding, this allows to consider 
N as a time variable, increasing with proper time t. We derive an equation 
for w from (15161) : 




The equations system dSHHl) is autonomous since w can always be written 
as a function of pi and 2 / 2 • This last statement needs explanations. It does 
not imply any coupling between matter, radiation and dark energy. Basi¬ 
cally, it is due to the fact that w can always be considered as a function 
of the redshift z. Then, since from (01151) . we get z = {p m o/ProU 2 /yi ~ 1), 
we can write, at least formally, w = w ( 2 / 2 / 2 / 1 ) • Finding the explicit form of 
w(yi,y 2 ) is obvious when the equation of state w is given as a function of 
z. For instance the linear equation of state|6] is w = wo + w\ z, and thus 
w(yi,y 2 ) = w{y 2 /yi) = w 0 + w 1 (p m0 /p r oy 2 /yi - 1). But w can also take 
more general forms than that of a function of y 2 /yi . For instance, this is the 
case when one looks for w mimicking the Dvali-Gabadadze-PorratifT] theory. 
Its modified Friedmann equation is H 2 ± H/r c = k/(6r c )(p m + p r ) with r c 
the crossover scalefS]. If we put the ±H/r c term on the right-hand side of 
the Einstein equations, defining the constant K by k = 2 r c K, this theory 
is mimicked by a dark energy theory with p = +2(r c K)~ 1 H (we take the 
”+” branch of the DGP theory since we assume p > 0). Using equation (HI), 
one shows that the corresponding w is then w = {y 2 — 3)/ [3(1 + 2/1 + 2 / 2 )] = 
w ( 2 / 1 , 2/2)0- Hence, in all this paper, we will write w as w(y\,y 2 ). 

Summarising, we assume that Universe is expanding and dark energy den¬ 
sity is positive. The field equations (1HU1) reduce to a system of two equations 

1 Even for the DGP theory, it is always possible to consider w as a function of the 

redshift and to write formally w = 10(2/2/3/1). In this case the mathematical relation 
10(3/2/1/1) = 10(2/1,2/2) expresses the fact that the two independent functions 2/1 and 1/2 
both depend on the redshift z. Indeed, since 1/1 = 2/i(z) and 1/2 = 1/2(2), it is always 
possible to write formally 1/1 (2/2) or equivalently the relation 10(3/2/2/1) = 10(3/1,3/2). 
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for y i and y 2 with u> (2/1,2/2) and can be studied with usual dynamical sys¬ 
tem method. Stability of the equilibrium points depends on the signs of the 
eigenvalues 



with 


A = -1 + 32/2 + 3(2 - 32 /i - 3 / j 2 )w ~ 3 (-l + 2/1 + 2/2X2/21%. + y\w yi ) 



and w Vi , the derivative of 10(2/1,2/2) with respect to 2/2- 

2.1 Equilibrium points 

Now, we look for the equilibrium points of the system ([SK]) whatever w. 
First we show that there is no equilibrium point with both 2/1 7^ 0 and 
2/2 / 0 . Let us assume the opposite, that is there is an equilibrium point in 
2/1 = 2/10 + 0 and 2/2 = 2/20 / 0 . Then, y\ will vanish if w = 3( _ 1+ y° 0+;;2o) ■ 
However, introducing this last expression in ([HD, it implies that y 2 will vanish 
only if 2/20 = 0 . This contradicts our assumption that 2/2 at equilibrium 
should not be zero. The same thing arises if we determine w from y 2 and 
try to make vanish y\ since then 7/10 should be equal to zero at equilibrium. 
Consequently, the equilibrium points are all such that 2/1 = 0 or/and 2/2 = 0 . 
Firstly, we consider the case y 1 = 0 . Then y\ disappears if w is such that 
2/i(l — 2/1 — 2/2XI + w) = 0 . This is always the case when w is finite. But if 
w diverges, this equality have to be true otherwise there is no equilibrium 
point on the line 2/1 = 0 . y 2 disappears in the following cases, taking into 
account that 0 < 2/2 < 1 

•2/2 = 0 with w such that 2/2(1 — 2/1 — 2/2 )(1 + w) = 0 

• 2/2 = 2/20 / 0 or 1 with w = 1/3 

• 2/2 = 1 with w such that (1 — 2/1 — 2/2) (1 + w)= 0 

Secondly, we consider the case 2/2 = 0 . Then y 2 disappears if w is such that 
2/2(1 — 2/1 — 2/2XI + w) = 0 . 2/1 disappears in the following cases, taking into 
account that 0 < 2/1 < 1 

• 2/1=0 with w such that 2/1 (1 — 2/1 — 2/2) (1 + w)= 0 

• 2/1 = 2 /io / 0 or 1 with w = 0 

• 2/1 = 1 with u; such that (1 — 2/1 — 2/2) (1 + w) = 0 
Summarising, we thus have the following equilibrium points 
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• (j/i) 2/2) = ( 1 , 0 ): this is the usual matter dominated equilibrium point. 
We call this point M. When w diverges, there is a M point if (1 — 

2/1 ~ 2/2 )w does not diverge. 

• (2/1! 2/2) = ( 0 , 1 ): this is the usual radiation dominated equilibrium 
point. We call this point R. When w diverges, there is a R point if 
(1 — ?/i — 2/2)10 does not diverge. 

• (2/1,2/2) = ( 0 , 0 ): this is the dark energy dominated equilibrium point 
for which Universe expansion can be accelerated (see below). We call 
this point D. When w diverges, there is a D point if y\w and y^w do 
not diverge. 

• (2/112/2) = (2/101 0 ) with i/io a non vanishing constant such that w(yio, 0 ) = 
0 : at this equilibrium point, dark energy mimics cold dark matter and 
we call it DM. Around DM , where w ~ 0 , dark energy mimics warm 
dark matter[T 7 j. In the rest of the paper we will not distinguish be¬ 
tween these two cases and we will write that in (when DM is a source 
or a sink for a trajectory in the phase space) or near (when DM is a 
saddle for a trajectory) DM, dark energy mimics matter and radiation 
is negligible. If there is no value 2/10 such that w(yio, 0 ) = 0 , there is 
no equilibrium point of this kind. Otherwise there are as many equi¬ 
librium points of this kind as there are constant values of y\ making 
w(yi, 0 ) vanishing. 

• (2/i)2/2) = ( 0 , i/20) with 2/20 a non vanishing constant such that w( 0 ,2/20) = 
1 / 3 : at this equilibrium point, dark energy mimics radiation and we 
call this kind of point DR (as for DM, in the rest of the paper we will 
not distinguish the cases when DR is a source or a sink for a trajectory 
with w = 1/3 or when it is a saddle and w ~ 1/3 during a transient 
period of time). If there is no value 2/20 such that w( 0 , 2/20) = 1 / 3 , there 
is no equilibrium point of this kind. Otherwise there are as many equi¬ 
librium points of this kind as there are constant values of 2/2 making 
MO,2/2) equal to 1 / 3 . 

2.2 General considerations about equilibrium points stabil¬ 
ity 

It is not possible to fully study the stability of the above equilibrium points 
without specifying w. However we can get some general and useful results. 

• The M point cannot be a source since we assume Universe is expanding 
and thus at early time, radiation always dominates matter. However, 
it can be a sink or a saddle. For instance when w and its derivatives 
w yi and w V2 do not diverge in M, (| 8 j) indicates that it is a saddle point 
when rc(l, 0 ) < 0 and a sink when re(l, 0 ) > 0 . 
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• The R point cannot be a sink since we assume Universe is expanding 
and thus at late time, matter always dominates radiation. However, it 
can be a source or a saddle. For instance, when w and its derivatives do 
not diverge in R, (JHJ) indicates that it is a saddle point when w(0, 1) > 
1/3 and a source when u>(0,1) < 1/3. 

• The D point can be a saddle, a source or a sink. For instance, when w 
and its derivatives do not diverge in D, ([8]) indicates that it is a saddle 
point when 0 < w( 0, 0) < 1/3, a source when w( 0, 0) > 1/3 and a sink 
when w( 0, 0) < 0. 


• Concerning the DM points, since we assume an expanding Universe, 
they can never be a source because in DM, 2/2 / 2/1 oc 1 + z tends to 
zero. However they can be a sink or a saddle. For instance if w' s 
derivatives do not diverge in DM, then from (O it comes that DM 
are saddle points when J^(?/io, 0) > 0 and sinks when ^(yiO)0) < 0. 


• Concerning the DR points, since we assume an expanding Universe 
, they can never be a sink because in DR, 2/2/ 2/1 oc 1 + z diverges. 
However they can be a source or a saddle. Indeed if re’s derivatives 
do not diverge in DR, using ([HI), we deduce that DR are saddle points 
when |^(0,y 20 ) < 0 and sources when ^(0,y 2 o) > 0. 

There is no limit to the number of equilibrium points of the system B 
0. For instance, if w = Ui=i(yi ~ Vu) YYjLi(y2 ~ 2/2 j) + 1/ YIi=x(-yii)/3 , 
there are R, M and D equilibrium points but also n equilibrium points DMi 
in (y 1 , 1 / 2 ) = (yn,0) with w = 0, yu some non vanishing constants and m 
equilibrium points DRj in (yi,y 2 ) = (0,y 2j ) with w = 1 / 3 , y 2 j some non 
vanishing constants. The case n = 3 and m = 1 is plotted on figure [T[ 

From a physical viewpoint, all source (sink) points are such that y\ = 0 
(respectively y 2 = 0). Thus, Universe always starts without matter and ends 
without radiation, even mimicked by dark energy. Moreover, by definition, 
Universe expansion accelerates when d 2 a/dt 2 = dH/dt + H 2 > 0, that is 


1/2 [-3 - y 2 + 3w(—l +y\+ y 2 )\ + 1 > 0 


Without surprise this last expression shows that the equilibrium points M 
and R cannot generate an accelerated expansion whatever the behaviour of 
w (which cannot diverge faster than (—1 + y\ +y 2 )~ 1 near M and R). DM 
and DR cannot generate an accelerated expansion neither. This is finally 
only possible for the equilibrium point D for which (yi,y 2 ) = (0,0) and 
when w < —1/3, a well known result 


2 Other equilibrium points with accelerated expansion could exist with a negative energy 
density for dark energy but we do not consider this case in that paper 
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Figure 1: Phase plane for w = (yi — l/4)(yi — l/2)(yi — 3/4) [( 2/2 — 1/3) — 32/9]. There are 
three equilibrium points in y% = 0 such that w = 0 with y\ = 1/4, 1/2 and 1/3 and one equilibrium 
point on y\ = 0 such that w = 1/3 with y 2 = 1/3. 

Last, we examine the dominant energy condition at each equilibrium 
point. For dark energy, it writes | p \ /p =| w |< 1. This condition can be 
violated in M, R and D, in particular in the case of a ghost dark energy and 
when w diverges. The condition always holds in DM and DR since then 
w is 0 and 1/3 respectively. If now we consider matter, radiation and dark 
energy as a unique fluid of density ptotai = p + p m + Pr, then the dominant 
energy condition writes | p tota i | / p to tai =| l/3y 2 +w(l-yi-y 2 ) |< 1 or ^ = 
-k/2(p tota i+p tota i) = 1/2 H 2 [—37/1 - 4 y 2 - 3(1 + w)(l-yi~ y 2 )\ < 0 when 
Ptotai < 0. These inequalities are always true at the points M, R, DM and 
DR, even when w diverges in M and R since then equilibrium is reached 
only if w « 1/(1 —y\ — 2 / 2 ) • The dominant energy condition is only violated 
in D, once again when | w |> 1. To be more specific about ^/H 2 , we note 
that its values in R, DR, M, DM and D, are respectively —2, —2, —3/2, 
—3/2 and —3/2(1 + w). The acceleration parameter q = 1 + has 

thus the same negative values in R and DR or in M and DM, in agreement 
with the fact that these two couples of equilibrium points have respectively 
the dynamics of some radiation and matter dominated epochs. 

3 Universe dynamics 

In this section, we look for some sequences of epochs describing a Universe 
starting by some early radiation dominated epochs, followed by some matter 
dominated epochs and ending by an accelerated expansion. Such sequences 
are described in the phase space ( 2 / 1 , 2 / 2 ) by trajectories having a source in 
R or DR, passing near some saddle points R or DR, then near some saddle 
points M or DM and ending in D. For instance, a phase space trajectory 
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starting in R, repelled by a first saddle point DR, then a second saddle 
point DM and ending in D will be noted R —>■ DR —>■ DM —>• D. There 
are an infinite number of possibilities depending on the form of w. However 
we will show in the next section that if w is not singular at some saddle 
equilibrium points, there are only ten of such sequences. We present them 
in this section by using polynomial forms of w(y 1 , 2 / 2 )- They are such that 
w and its derivatives stay finite whatever 0 < jji < 1. This also allows to 
use results of section 12.21 and to get easily some w reproducing the above 
desired sequences of epochs. 

3.1 Three successive epochs 

In this subsection, we consider the above described sequences with three suc¬ 
cessive epochs. Without surprise, there are four possibilities that we name 
£31 to £34 (where in Sij, i stands for the number of successive epochs in a 
sequence, and j is the sequence number): 

• S31: R^M ^D 

This is the sequence followed by the AC DM model illustrated on figure [2] 
and most of the dark energy models described in the literature. 



0.0 0.2 


.8 1.0 


Figure 2: Phase plane for the cosmological constant w = — 1. Grey area is the zone of the 
phase plane where Universe expansion is accelerated. The points are the equilibrium points. A 
trajectory is plotted as a dotted line. 


• £ 32 : R -> DM -A D 

Dark energy mimics dark matter at intermediate time and accelerates ex¬ 
pansion at late time. There is thus at least one value yio of y\ such that 
w(y 10 , 0) = 0. Although one could think that something similar occurs with 
the Chaplygin gas)10j. this last theory rather follows the £31 sequence with 
w = 0 in ( 2 / 1 , 2 / 2 ) = (0,1), i.e. at early times. An example of sequence £32 is 
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shown on figure 03 



Figure 3; Phase plane for w = —1 + 1.42yi with a saddle equilibrium point DM in fi/ 1 , 2 / 2 ) = 
(0.7,0). Grey area is the zone of the phase plane where Universe expansion is accelerated. The 
points are the equilibrium points. A trajectory is plotted as a dotted line. 


• S 33 : DR^M^D 

Dark energy mimics radiation at early time and accelerates expansion at late 
time. There is thus at least one value 2/20 °f V 2 such that w(0,y2o) = 1/3. 
This kind of dark energy plays the same role as Chaplygin gas at early time 
but mimicking radiation (w —>• 1/3) instead of matter. An example of such 
a sequence is shown on figure [U 



0.0 02 


Figure 41 Phase plane for w = — 1 + 7.98^2 — 6.58^2 with a source equilibrium point DR 
in ( 2 / 133 / 2 ) — (0,0.2). Grey area is the zone of the phase plane where Universe expansion is 
accelerated. The points are the equilibrium points. A trajectory is plotted as a dotted line. 
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• 5 34 : DR -> DM -> T> 

Universe evolution is driven by dark energy that mimics radiation at early 
time, (dark) matter at intermediate time and accelerates expansion at late 
time. In this case, there are at least one value 2/10 of y\ and one value 7/20 of 
2/2 such that respectively 7 / 7 ( 2 / 10 ,0) = 0 and 7 / 7 ( 0 , 2 / 20 ) = 1/3- An example of 
such a sequence is shown on figure [5j 



Figure 5: Phase plane for w = —1 + 1.25yi + 1.90 j/ 2 with a saddle equilibrium point DM in 
(y 1,7/2) — (0.8,0) and a source equilibrium point DR in (yi,y 2 ) = (0,0.7). Grey area is the zone 
of the phase plane where Universe expansion is accelerated. The points are the equilibrium points. 
A trajectory is plotted as a dotted line. 


To our knowledge, cosmological models following sequences S 32 to S34 have 
not been studied in the literature. Indeed, most of dark energy models fo¬ 
cused on late time epoch. Consequently the role plays by a dark energy as 
an additional source of radiation at early times or/and of dark matter at 
intermediate times, is usually not envisaged. 

3.2 Four successive epochs 

In this subsection, we consider sequences with four successive epochs. There 
should be 12 possibilities but only 6 of them exist when w is not singular at 
some saddle equilibrium points (see section |4|). 

• S 41 : DR^ M -> D 

This first possibility is similar to usual dark energy models but there are 
now two radiation epochs at early times. For the first one, there must be 
a value 7/20 of 2/2 such that w(0, 2/20 ) = 1/3 and then, dark energy mimics 
radiation. An example of such a sequence is plotted on figure [6j 

• S 42 : R-* DR^ M -> D 
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Figure 6 ; Phase plane for w = —1 + 3.33y2 with a source equilibrium point DR in (yi,y 2 ) = 
(0,0.4). Grey area is the zone of the phase plane where Universe expansion is accelerated. The 
points are the equilibrium points. A trajectory is plotted as a dotted line. 


This sequence is similar to S 41 but the two radiation epochs are reversed. 
An example of such a sequence is plotted on the first panel of figure [71 



yi yi 


Figure 7; Phase plane for w = — 1 + I0y2 — 16.66y2 with two equilibrium points DR , one saddle 
in (yi, y 2 ) = (0, 0.4) and one source in (yi, y 2 ) = (0, 0.2). Grey area is the zone of the phase plane 
where Universe expansion is accelerated. The points are the equilibrium points. Two trajectories 
are plotted as dotted lines. The first panel shows a R —> DR —> M —> D trajectory and the second 
panel a DR —> DR —>• M —> D trajectory. 


• S 43 : DR-> DR^ M -> D 

During the two radiation epochs dark energy mimics radiation but in dif¬ 
ferent proportions. This time we thus need two non vanishing values of 7/2 
such that w( 0 , 7 / 2 ) = 1/3. An example of such a sequence is plotted on the 
second panel of figure |71 
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• S44: DR — > R —» DM —>■ Z) 

Dark energy mimics radiation during the first radiation epoch. The second 
epoch is dominated by pure radiation. During the third epoch, dark energy 
mimics matter. Thus there should be at least one value yio of y\ and one 
value 2/20 of y 2 such that respectively w(y 10 , 0) = 0 and u>(0, 2 / 20 ) = 1/3. An 
example of such a sequence is plotted on figure 0 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 



Figure 8; Phase plane for w = — 1 + 0. 62yi +0.78?/^ — 5.312/2 +10.3I2/2 with a saddle equilibrium 
points DM in (271,2/2) = (0.8,0) and a source equilibrium point DR in (271,2/2) — (0,0.7). Grey 
area is the zone of the phase plane where Universe expansion is accelerated. The points are the 
equilibrium points. A trajectory is plotted as a dotted line. 


• S 45 : R DR DM -A D 

This is similar to S44 but with the two radiation epochs reversed. An ex- 
ample of such a sequence is plotted on figure [9l 

• S 46 : DR —>■ DR —»■ DM -»• D 

Dark energy mimics radiation during the first and second epoch, matter dur¬ 
ing the third epoch and finally accelerates expansion at late time. All Uni¬ 
verse dynamics is thus driven by dark energy. At least two values of j /2 have 
to be such that w(0,y2) = 1/3 and one value of y\ such that w(yi,0) = 0. 
An example of such a sequence is plotted on figure [TUJ 

In the next section, we show that these 10 sequences, that describe a Uni¬ 
verse starting with some radiation epochs followed by some matter epochs 
and ending by an accelerated expansion, are the only ones that can be got 
if w is not singular at some saddle equilibrium points. 
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Figure 9; Phase plane for w = —3/2 — 3yi + 17.66y^ — 14.1+ 6 . 33 y 2 — 5 . 33^2 with two 
equilibrium points DM, one sink in ( 2 / 1 , 3 / 2 ) — (0.85,0) and one saddle in ( 2 / 1 , 3 / 2 ) = (0.6,0), and 
two equilibrium points DR, one source in ( 3 / 1 , 3 / 2 ) = (0,0.5) and one saddle in ( 3 / 1 , 3 / 2 ) = (0,0.68). 
Grey area is the zone of the phase plane where Universe expansion is accelerated. The points are 
the equilibrium points. A trajectory is plotted as a dotted line. 



yi 


Figure 10: Phase plane for w = -1.12-5.74yi+23.12y2_i7.29yf+ 11.61y 2 -24.70y| + 15.72^ 
with two equilibrium points DM, one sink in ( 3 / 1 , 3 / 2 ) = (0.86,0) and one saddle in ( 3 / 1 , 3 / 2 ) = 
(0.59, 0), and three equilibrium points DR, two sources in (y\, 7 / 2 ) = (0, 0.76) and (yi, 2 / 2 ) — (0, 0.2) 
and one saddle in ( 3 / 1 , 3 / 2 ) — (0,0.6). Grey area is the zone of the phase plane where Universe 
expansion is accelerated. The points are the equilibrium points. A trajectory is plotted as a 
dotted line. 


4 Discussion 

In the previous section, we reproduced 10 sequences with no more than four 
epochs describing Universe starting with some radiation epoch(s) followed 
by some matter epoch(s) and ending by an accelerated expansion. But, of 
course, there are an infinity of them, with as many epochs as desired. In¬ 
deed one can always choose any form of the scale factor a and then derive 
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the corresponding w. We remark that a common property of this infinity 
of sequences, but the ten we found in the previous section, is that they al¬ 
ways contain two successive saddle epochs dominated by matter (noted as 
DM -H- M or DM -H- DM) or radiation (noted as DR -H- R or DR -H- DR). 
In the same time we show in the appendix that when w and/or its derivatives 
with respect to y\ and 1/2 are finite, it is not possible to get two successive 
saddle epochs dominated by matter or radiation. Consequently, we deduce 
that to get additional sequences of epochs with respect to the ten we found 
in section O w and/or its first derivatives with respect to y\ and/or y2 have 
to diverge at some saddle equilibrium points of type M, DM, R or DR and 
thus be singular at some of these points. From a physical viewpoint, since a 
trajectory pass near a saddle point and not exactly through it, it means that 
the w corresponding to this trajectory will have some pronounced bumps or 
dips. 

Such a behaviour seems not favoured by observations. Hence in m , no 
evidence is found for a spike in the equation of state at early time that 
would have left an imprint on the CMB anisotropy pattern and the rate of 
growth of large scale structure. In [14] it is shown that large fluctuations of 
the equation of state are allowed only at ultra-low redshifts z < 0.02. An¬ 
other way to interpret the divergence of w and its derivatives is to consider 
the adiabatic sound speed. It writes 

9 dp/dt w' 

ct = ——— = w -7-r 

dp/dt 3(1 + w) 

or with the normalised variables 

2_ V2(w y2 - w yi yi - w m y 2 ) + 3 w(-l + 2/1 + y 2 ){w yi y l + w m y 2 ) 

C °~ W+ 3(1 + w) 

When Cg > 1, perturbations move at superluminal speeds[T5] whereas when 
(? s < 0 , the system is perturbatively unstable on small scales and therefore 
unphysical. Then, as explained in m, the adiabatic pressure would accel¬ 
erate the collapse of density perturbations unless entropic perturbations be 
taken into account and the the sound speed of perturbations be no more the 
adiabatic sound speed. Also, we are going to examine the values of at 
each equilibrium point to check when 0 < < 1 . 

Firstly, we assume that w and its derivatives w Vi do not diverge. Then 

= w for all the equilibrium points. It will be negative in M and D when 
the first point is a saddle and the second one a sink since then w < 0 (see 
subsection 12.21) . In R, there is no restriction on the sign of w when this 
point is a source or a saddle. However, we may have < 0 only when R is 
a source or c? s > 1 only when R is a saddle. There is no problem for in 
DM and DR since w = 0 and w = 1/3 respectively. 
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Secondly, we assume that w does not diverge but its derivatives w Vi . Then 
| c 2 |>> 1 in M when the w yi diverge faster than l/y 2 , in DM when 
w Vi » l/y 2 and w yi » 1/w, in R or DR when the w yi diverge faster 
than I/ 7/1 or finally in D when w yi diverges faster than 1 /t ji with i = 1,2. 
If these inequalities are not respected, then the divergence of the w yi do not 
intervene in the value of | c 2 | and we recover the previous results, when the 
w yi do not diverge. 

Last, when w diverges at equilibrium, generally c 2 always diverges (but in 
the special case for which the divergence of w is compensated by this of the 
second term in the right side member of equation (|9|)). 

Summarising, if we assume that M is a saddle and D is a sink, c 2 is always 
negative at these points, a well known fact for D, underlined in ESJ. For 
the other equilibrium points, we will have 0 < c 2 < 1 only if w does not 
diverge and then in R if w yi « 1/yi and 0 < w < 1, in DM if w yi « l/y 2 
and w yi « 1/w, in DR if w yi « ljy\. Hence divergence of w and/or 
its derivatives w yi generally reduces the number of equilibrium points where 
the adiabatic sound speed is well defined whereas we will always have to 
take into account entropic perturbations and to consider an effective sound 
speed [16] at intermediate and late times in M and D respectively, to study 
structures formation in presence of dark energy. 

To illustrate the fact that a form of w with singular behaviours at some 
equilibrium points M, R, DM and DR allows to get some sequences different 
from the ten we present in the previous sections, we try many forms of w with 
such properties. This is generally not enough to get these new sequences. 
This implies that a singular w at some equilibrium points is a necessary but 
not sufficient conditions to get them. Finally, we found that the function w 
given by 

w = -(i/i - 0.3) 2 (t/i - 0.7) 2 /(( 7/1 - 0.3) 2 + t/ 2 )/((t/i - 0.7) 2 + t/ 2 ) - y 2 /yi 

allows to reproduce a new sequence, the sequence R —>• DM —>• DM —>• D. 
w is zero in (t/i, 7 / 2 ) = (0.3,0) and ( 7 / 1 , 7 / 2 ) = (0.7,0) but its derivatives are di¬ 
verging at these points (see below). The corresponding phase space is plotted 
on figure [TTJ A sequence R —>• DM —>• DM —>• D is shown as a dotted line 
but is difficult to observe clearly. Also we plotted on figure [T2] the evolution 
of 7/1 (IV) and w(N) when 7/1 (0) = 0.27 and 7 / 2 ( 0 ) = 8.27 x 10 -5 , the present 
day observational values for matter and radiation density parameters, w has 
two pronounced (but not singular since the trajectory pass close but not ex¬ 
actly on the DM points) bumps as shown by the second graph on figure [T2l 
These bumps correspond to the trajectory approaching the two DM saddle 
points in ( 7 / 1 , 7 / 2 ) = (0.3,0) and (t/i, 7 / 2 ) = (0.7,0), where w is zero and its 
derivatives are diverging. If these bumps are so pronounced, this is due to 
the fact that radiation density parameter should be very small during matter 
era. Hence the trajectory pass very close to the points where w is singular. 
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Figure 11: Phase plane for w = —(j/i—0.3) 2 (yi—0.7) 2 /((yi —0.3) 2 )/((yi —0.7) 2 +j/§)—J/ 2 / 2/1 
with two DM equilibrium points in ( yi,y 2 ) = (0.3,0) and (f/ 1 , 2 / 2 ) = (0.7,0). The first one is a 
saddle. The second one is a saddle when y\ <0.7 and a sink otherwise. Grey area is the zone of 
the phase plane where Universe expansion is accelerated. The points are the equilibrium points. 
A trajectory presenting a sequence R —> DM —> DM —> D is plotted as a dotted line. 

Note that between the saddle points, the expansion undergoes a transient 
acceleration with w = — 1 that does not correspond to an equilibrium state. 



N N 

Figure 12: yi(N) and ro(JV) when w = — ( 3/1 - 0.3) 2 (yi -0.7) 2 /((m -0.3 ) 2 + 2 /§)/(( 2 /i -0.7 ) 2 + 
j/ 2 ) — 3 / 2 / 2/1 with j/i(0) = 0.27 and 3 / 2 ( 0 ) = 8.27 X 10 -5 . The two bumps in the evolution of w 
correspond to the trajectory in the phase space approaching the DM saddle points. 


To conclude, in this paper we looked for sequences of epochs able to repro¬ 
duce our Universe history when gravity is described by General Relativity 
with radiation, matter and a dark energy defined by an equation of state w. 
We assumed that Universe is expanding and dark energy density is positive. 
Using dynamical systems methods and a set of normalised variables, we de¬ 
fined all the possible equilibrium points of this theory with these variables 
and gave some general properties about their stability. These equilibrium 
points describe some classical matter or radiation dominated epochs but 
also dark energy epochs during which dark energy mimics matter, radiation 
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or accelerates Universe expansion. When w is not singular at some saddle 
equilibrium points, there are ten possible sequences of epochs describing 
our Universe history. There are four sequences with three epochs and six 
sequences with four epochs. In the dark energy literature, cosmological mod¬ 
els often follow sequence S 31 . Here we find some new possibilities for which 
dark energy can also mimic radiation in a similar way it mimics matter in 
the Chaplygin model. Although these possibilities may not be distinguished 
by observations like supernovae, they could be of interest from structure 
formation viewpoint. 

To get more than these ten sequences, it is necessary but not sufficient 
that w or/and its first derivatives with respect to y\ and/or y^ diverge at 
some equilibrium points M, DM, R or/and DR, i.e. at intermediate time. 
However such properties for w seem to be discarded by observations[13] [13j. 
These properties also increase the number of equilibrium points where the 
adiabatic sound speed is ill-defined. Moreover, the divergence of w leads 
to the violation of the dominant energy condition. All these clearly distin¬ 
guish the ten above sequences from an infinity of ways to describe Universe 
expansion. 


Appendix - Special case: w{y\, 7 / 2 ) and its derivatives 
are finite 

Here we show that when w and its derivatives are finite, we cannot have two 
successive saddle epochs dominated by matter or radiation, i.e we cannot 
have a transition between equilibrium points like DM -H- DM, DR -H- DR, 
DM -H- M or DR -H- R. For that we use the results of section 12.21 where the 
stability of the equilibrium points is given in such a case. 

Let us consider the equation of state 10 ( 2 / 1 , 2 / 2 )- 10 and its derivatives with 
respect to y\ and 2/2 are finite everywhere on the phase plane. We define 
the intersection of the surface w(yi, 2 / 2 ) with the plane 2/2 = 0 as a C l curve 
(Cm) containing all the equilibrium points with 2/2 = 0, i.e. the points D, 
M and n points of type DM. We call these last successive points DM\, 
DM2-..DM n . DM\ is the nearest equilibrium point from D and DM n the 
nearest equilibrium point from M. This is illustrated on figure [13] with 
n = 4. In the same way, we define (Cr) the intersection of the surface 
io(2/i; 2 / 2 ) with the plane y\ = 0. It thus contains all the equilibrium points 
D, R, DR\...DR m with y\ = 0. 

First, we show that when DM % is a saddle then DMi + \ is a sink and 
DMi- 1_2 is a saddle again. When DM % is a saddle, w = 0 and ^(yio,0) > 0 
at this point. Following the (Cm) curve, it will necessarily cross the w = 0 
value again in DMi + \ with this time jj^(yio, 0) < 0 like on figure [13l Hence 
DM{. 1-1 will be a sink. The same reasoning leads to show that in DMi + 2 , 
^-(2/io,0) > 0 and it will be a saddle. Hence we have shown that when 
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Figure 13: The intersection of a surface w(yi ■ y 2 ) with the plane 3/2 = 0 forms a curve (Cm) 
containing all the equilibrium points with 3/2 — 0. 

there are several DM points on r /2 = 0, if the first one is a saddle (sink), 
the next one is a sink (respectively saddle) and so on. In the same way one 
can prove that if there are several DR points on y\ = 0, if the first one is 
a saddle (source), the next one is a source (respectively saddle) and so on. 
It follows that when w and its derivatives are finite everywhere, we cannot 
have two successive saddle epochs like DM, —>• DMj because any trajectory 
going from one saddle point DM. t to the other saddle point DMj would 
isolate a sink point between DMi and DMj from any source. One shows in 
the same way that the transition between two successive saddle epochs like 
DRi —» DRj is not possible because any trajectory going from one saddle 
point to the other would isolate a source point between DRi and DRj from 
any sink. 

Now we want to show that if DM n is a saddle point, the next equilibrium 
point on a trajectory cannot be a saddle point of type M and vice-versa. 
Indeed if the point DM n is saddle, then in DM n , w = 0 and ^(y leu 0) > 0. 
Following the continuous curve (Cm), we thus arrive in M with w > 0 0 . 
Since we assume that w and its derivatives are finite, we know from section 
I2.2l that M is then a sink and cannot be a saddle. Consequently, a transition 
from a saddle point DM n to a saddle point M is impossible. In the same 
way one shows that a trajectory cannot go through a saddle point M to a 
saddle point DM n . 

Last we show that if DR n is a saddle point, the next equilibrium point 
cannot be a saddle point of type R and vice-versa. Indeed if DR n is a 
saddle, then in DR n , w = 1/3 and ^(0, ?/ 2 o) < 0. Following the continuous 
curve (Cr), we thus arrive in R with w < 1 / 30 . Once again, we derive from 
section [2T2l that R is then a source and never a saddle, w and its derivatives 
being finite. Hence, a transition from a saddle point DR n to a saddle point 
R is impossible. In the same way one shows that a trajectory cannot go 

3 Otherwise, w would be equal to zero between DM n and M on (Cm) and there will be 
another DM point closest to M than DM n which would contradict our initial assumptions 

4 Otherwise, w would be equal to zero between DR n and R on (Cr) and there will be 
another DR point closest to R than DR n which would contradict our initial assumptions 
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through a saddle point R to a saddle point DR n . 

We conclude that when w and its derivatives are finite, we cannot have two 
successive saddle epochs dominated by matter or radiation. 
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